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INTRODUCTION 


The  problem  we  wish  to  investigate  is  that  of  the  elastic  stability 
of  thin  circular  conical  shells  under  uniform  external  pressure.  Both 
conical  frustra  and  complete  cones  are  to  be  considered. 

In  those  cases  of  shell  stability  where  it  is  difficult  or  impossible 
to  solve  the  problem  exactly,  the  direct  methods  of  the  variational  calculus 
deserve  special  attention.  The  more  complex  the  .problem,  the  more  advan¬ 
tageous  is  the  use  of  such  methods. 

Several  theoretical  investigations  of  the  problem  of  conical  shell 
stability  employing  an  energy  procedure  have  been  carried  out  in  recent 
years  without  resolving  the  problem.  One  of  the  underlying  problems  of 
stability  analyses,  and  the  probable  source  of  much  of  the  disagreement, 
is  that  they  are  beset  with  considerable  algebraic  and  numerical  difficul¬ 
ties  which  necessitate  the  introduction  of  approximations.  None  of  these 
approximate  theories  yield  results  in  close  agreement  with  one  another  or 
with  experiment  over  a  wide  range  of  geometries. 

In  the  interest  of  simplicity  most  investigators  simplify  the  problem 
by  neglecting,  in  the  strain  energy  expression,  those  displacement  terms 
of  cubic  and  higher  order,  and  employing  the  principle  of  virtual  displace¬ 
ments.  This  states,  "Of  all  the  displacements  satisfying  given  boundary 
conditions,  those  which  also  satisfy  the  equilibrium  conditions  make  the 
total  potential  energy,  tt,  assume  a  stationary  value.  "  In  variational 
terms,  this  corresponds  to  setting  the  first  variation  of  the  total  poten¬ 
tial  energy  equal  to  zero. 

However,  such  a  procedure  minimizes  the  influence  of  rotation  terms 
on  the  deformation.  In  thin  flexible  bodies  such  as  shells,  this  implies 
an  effective  stiffening  of  the  shell  and  a  higher  estimate  of  the  critical 
load. 


The  object  of  this  investigation  is  to  present  a  systematic  procedure 
to  account  more  adequately  for  the  influence  of  rotation  effects  in  thin 
shells.  We  first  introduce  a  form  of  strain  displacement  relations  in 
which  specific  emphasis  is  placed  on  the  rotation  terms.  Next,  higher- 
ordered  terms  are  retained  in  the  expression  for  the  total  potential  energy 
of  the  system,  and  an  appropriate  stability  criterion  is  formulated. 

Now,  if  a  functional  such  as  tt  is  quadratic  in  its  independent  vari¬ 
ables  and  their  derivatives,  a  stationary  value  is  unique  provided  it 
exists.  (The  Euler  equations  and  the  natural  boundary  conditions  are 
linear. )  If  the  functional  is  of  higher  order  in  its  independent  variables, 
multiple  stationary  values  may  exist.  (The  Euler  equations  are  nonlinear.  ) 
This  introduces  the  possibility  of  alternative  equilibrium  positions  some 
of  which  are  stable,  others  not.  An  equilibrium  position  is  designated  as 
as  stable  when  the  virtual  work  is  negative  for  all  admissible  "virtual 
displacements,"  i.e.,  when  additional  external  forces  are  necessary  to 
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displace  the  system  from  its  equilibrium  position.  In  variational  nota¬ 
tion,  an  extremal  defining  a  stable  equilibrium  position  must  maximize 
the  work  functional  or  minimize  the  total  potential  energy  functional,  tt. 
Since  (for  small  virtual  displacements)  an  extremal  is  a  minimum  when  the 
second  variation  is  positive  definite  (i.e. ,  positive  for  all  admissible 
argument  functions  in  the  neighborhood  of  the  extremal),  the  criterion 
for  stable  equilibrium  is  8g77  >  o.  By  the  same  token,  an  unstable  equil¬ 
ibrium  position  is  characterized  by  <  o.  Hence  the  condition  for  the 
onset  of  instability  for  thin  sheila  is  =  o.  * 

The  theory  thus  formulated  is  valid  for  small  virtual  displacements 
about  the  equilibrium  position,  and  is  sometimes  referred  to  as  the  in¬ 
finitesimal  theory  of  buckling.  It  was  first  proposed  by  Trefftz*  and 
was  employed  in  the  study  of  the  buckling  of  cylindrical  shells  under 
external  pressure  by  Langhaar  and  Boresi.®  and  more  recently,  in  the 
analysis  of  the  buckling  of  a  cylinder  under  circumferential  band  loads 
by  Brush  and  Field. ^ 

The  procedure  outlined,  with  the  general  strain-displacements  given 
in  Appendix  A,  is,  then  generally  applicable  to  the  stability  of  thin 
shells.  Explicit  estimates  of  critical  loadings  may  be  made  with  the  aid 
of  the  Rayleigh-Ritz  method.  It  is  shown  in  the  body  of  the  report, 
however,  that  in  using  this  method  care  must  be  taken  to  select  a  set  of 
displacement  functions  which  are  realistic  over  the  entire  range  of  geom¬ 
etries  to  be  considered. 


DERIVATION  OF  POTENTIAL  ENERGY  OF  THE  SYSTEM 


We  consider  the  problem  of  the  elastic  stability  of  a  conical  shell 
to  be  one  in  the  nonlinear  theory  of  elasticity,  with  small  strains  and 
small  but  finite  deflections. 

An  appropriate  formulation  of  the  strain-displacement  relations  for 
such  a  problem  is  derived  in  Appendix  A  in  terms  of  the  thin  shell  coor¬ 
dinate  system  (a-^ya^yz)  ^  where  a-^ya,^  are  lines  of  principal  curvature  on 

the  shell’s  middle  surface,  having  principal  radii  of  curvature  R-j^  and  Rg, 
and  Lame  coefficients  H^^jHg  of  the  form 

Hi  >  Aid  •  VSi) 

•  Ajd  .  VRs)  . 

and  z  is  the  thickness  coordinate,  normal  to  (aj^,ag).  In  terms  of  the 
middle  surface  displacements  (u,v,w)  in  the  {a^,a,g,  z)  directions,  the 
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strain  displacement  relations  may  be  written  as: 


where 
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and,  finally, 


0  =  -  e^rj  i//  =  -  egg  X  =  ©ii  +  ©22 


We  may  further  simplify  these  relations  and  still  retain  the  same 
degree  of  approximation  already  inherent  in  them  by  observing  that  wg, 
the  component  of  rotation  about  the  z  axis,  will  be  very  small  compared 
to  the  remaining  components  of  rotation  and  wg.  This  is  so  because 
the  shell  is  very  stiff  ("massive")  in  its  own  plane.  Accordingly  we  may 
omit  the  nonlinear  terms  containing  ojg  in  the  above  expressions  for  the 
strains. 


We  may  specialize  to  conical  geometry  with  the  aid  of  Figures  lA 
and  IB.  For  conical  coordinates: 


tti 


=  f 


flt  g  ®  d 


h  = 


Ag  =  r  =  Rg  sin  a  *  f  cos  a 


(4) 


and  the  appropriate  nonlinear  strain-displacement  relations  are: 
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It  should  be  pointed  out  that  the  derivation  of  equations  1  involves 
no  order  of  magnitude  estimates  on  the  relative  values  of  the  displace¬ 
ments  themselves,  but  only  on  the  over-all  character  of  the  defor^tion. 
For  this  reason  it  is  believed  they  are  more  generally  applicable  to  the 
stability  problem  of  shells.  For  example,  in  formulating  Donnell's 
equation,  the  number  of  circumferential  waves  in  the  buckling  pattern  is 
assumed  to  be  large.  On  this  basis,  w  is  assumed  to  be  small  compared 


with 


902 


9w 

and  V  is  assumed  negligible  compared  with  ~ 

30 


The  strain  energy  of  a  conical  shell  is  given  by 

U  =  -  Ub 

where 


(6) 


TT  Sh  Btt  a  ( 

^  ”  S(  l-v2)  •()  o 


2  +  +  2v 


* 


(¥M 


rd0d^ 


Eh^  2w  J. 

”  84(l-v^)  4  4  *  2{1-v)K^/)  rd0d^ 


}  (7) 


where  E  =  Young’s  modulus,  v  »  Poisson’s  ratio  and  h  ■  shell  thickness. 
Finally,  the  potential  energy  of  the  system  is  given  by 


“  Urn  ^  Ub  -  Ue^t 


(8) 


where  change  in  potential  of  the  external  forces.  In 

i^pendix  B  we  have  shown  that  for  an  external  pressure  p. 


^ext  “  P  «■  //(^  +  u){v^(w^-v)-(R  +  vi)^(vg  +  w  +  r)}d^d0 


+  P 


sin  a  < 


21Ai 

3 


+ 


e(A2^  +  Ag) 


r  (9) 


where  A-^  and  Ag  are  the  cross-sectional  areas  at  the  top  and  bottom  of 
the  cone,  and  e  is  the  average  extension  of  the  ends  of  the  cone. 
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DETERMINATION  OF  DISPLACEMENT  FUNCTIONS 


The  Rayleigh-Ritz  procedure  is  a  convenient  technique  for  approxi¬ 
mating  the  solution  of  buckling  problems  where  exact  solution  is  difficult. 
It  results  in  an  upper  bound  estimate  for  the  critical  load.  The  method 
involves  assuming  a  form  for  displacements  which  satisfies  the  boundary 
conditions  and  which  contains  a  set  of  arbitrary  parameters.  These  para¬ 
meters  are  then  chosen  to  optimize  the  quantity  desired. 


Ideally,  the  chosen  set  of  displacement  functions  should  also  present 
a  reasonable  description  of  the  actual  displacements  occurring  during  the 
deformation  of  the  shell.  The  estimate  of  the  critical  load  is  fairly 
sensitive  to  the  choice  of  displacement  functions. 


Another  way  of  expressing  this  is  that  an  optimal  set  of  displacement 
functions  used  in  a  Rayleigh-Ritz  procedure  is  one  which  yields  a  good 
approximation  to  those  displacements  which  actually  satisfy  the  equilibrium 
equations  of  the  shell  and  yet  remain  simple  enough  in  form  to  be  of  prac¬ 
tical  benefit. 


Niordson®  and  Radkowski®* ’  employed  the  (cylindrical)  displacement 
functions 


Uj^  =  A  cos  n  6  cos 

v^  ®  B  sin  n  d  sin  / 

Wj^  -  C  cos  n  6  sin  ■i 


(10) 


which  satisfy  conditions  of  simple  support  at  both  ends  of  a  truncated 
cone  and  are  also  the  exact  solutions  for  a  cylindrical  shell  of  finite 
length.  They  apparently  yield  satisfactory  results  for  very  short  conical 
frustra.  In  terms  of  the  important  "taper  ratio”  parameter  yS  *  (1  cosa)/ro, 
these  functions  are  good  estimates  of  the  buckling  pattern  for  small 
compared  with  unity.  Niordson  limited  the  application  of  his  results  to 
/8  <  1/3.  Radkowski  retained  terms  neglected  by  Niordson  and  thereby  ex¬ 
tended  the  range  of  application  of  his  results,  but  found  that  as  /S 
approached  unity,  these  displacements  introduce  a  singularity  into  the 
expression  for  p^p.  Radkowski  used  a  physical  argument  to  extrapolate 
his  results  to  complete  cones.  Grigolyuk®  selected  the  displacement 
functions 


Un 


Wn 


Ar^  cosf— ]sin  n  6 

\T0) 


Br^  sin| — jcos  n  6 
Cr^  sin^— jsin  n  6 


(11) 
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which  satisfy  the  apex  conditions  for  a  complete  cone: 

u  =  v  =  w  =  u,  j.  =  v,  p»w,  j.  =  o  atr  =  o. 

However,  they  are  not  satisfactory  for  a  truncated  cone  unless  modified 
somewhat.  For  example,  a  more  satisfactory  form  for  Wj^  for  truncated  and 
complete  cones  might  be: 

P  (  r  -  ro\ 

W]^  =  Cr^  sin  ttI  sin  n  6 . 

The  buckling  displacement  functions  employed  in  the  present  investigations 
are  of  the  form: 


r  /r\2  1 

•V 

cos  (ar  +  b)  cos  n  6 

r  /  r  \2  “1 

sin  (ar  +  b)  sin  n  6  f 

sin  (ar  +  b)  cos  n  6 

\ _ 7  \  ^0  / 

> 

(12) 


where 


TT  TITO 

a  ®  -  b  - 


r-j^  -  ro  r-j^  -  ro 


and 


while: 


fl(^)  -  1,  °  as  ^  -  1 


fl(^)  -  o,  f2(^)  -  1  as  /8  -  o. 


It  may  easily  be  seen  that  this  choice  of  displacements  satisfies 
suitable  edge  conditions  for  both  complete  and  truncated  cones.  Moreover, 
for  suitable  f^  and  fg  it  should  reproduce  the  dependence  of  the  buckling 
pattern  on  the  slant  length  of  the  cone  -  viz. :  for  short  frustra  the 
dimpled  region  extends  over  a  much  larger  percentage  of  the  slant  length 
of  the  cone. 
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In  particular,  Figure  2  contains  a  plot  of  wi/wimax.  ^  function 

of  the  slant  length,  ^/-l,  for  several  values  of  /3.  In  this  graph  we  have 
selected 

=  7“^  >  ^2 

1  +  P 


Examination  of  Figure  2  points  out  why  one  would  not  expect  the  Niordson- 
Radkowski  displacements  (yS  =  o)  to  be  suitable  for  complete  cones  yS  “  1,  and 
conversely,  why  Grigolyuk's  functions  (yS  =  l)  would  not  be  suitable  as  yS-'o. 

On  the  other  hand,  Seide®  who  has  employed  a  power  series  solution 
of  the  equilibrium  equations  has  also  computed  the  displacement  w  which 
emerges  from  his  solution.  A  comparison  of  Figure  2  with  his  results 
indicates  very  good  qualitative  agreement. 


Our  displacement  functions  in  final  form,  then,  are: 
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.  (14) 

-  /3)^ 

sin  (ar  +  b)  cos  n  6 

> 

where  we  recall  that  Uq,  Vq  and  w^  are  to  be  determined  from  the  linear 
(pre-buckling)  theory. 

If  we  assume  that  the  pre-buckling  state  of  the  shell  is  adequately 
described  (away  from  the  edges)  by  the  membrane  theory  of  thin  shells, 
then  we  may  write  for  the  direct  stress  resultants  in  the  cone: 


N. 


pr 


,  N 


pr 


e 


^  2  sin  a  “  sin  a 

From  the  linear  strain-displacement  and  stress-strain  laws  we  have 


(15) 


Bu  1 
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so; 

(l  -  2v)  pr^ 

n  =  -  . 

4Eh  sin  a  cos  a 


(16) 


To  find  Wq  we  have: 


£/,  =  —  (Na  -  I'Ni; )  =  —  (w  sin  a 
^  Eh  ^  r 


u  cos  a)  , 


so: 


3pr^ 


4Eh  sin*  a 


(17) 


Due  to  rotational  symmetry  prior  to  buckling,  =  o. 


DETERMINATION  OF  THE  BUCKLING  LOAD 

The  stability  of  the  system  can  be  investigated  by  considering  a 
virtual  displacement  from  a  loaded  but  unbuckled  equilibrium  configuration 
and  examining  the  corresponding  second  variation  of  77. 

If  to  Uq,  Vq  and  Wq  (the  pre-buckling  equilibrium  configuration)  are 
added  the  (virtual)  displacements  u^f  v^  and  then  77  becomes: 


77  +  A77  =  77  + 


S  2i77  6  jj77  B  ^77 

"if  *  "sT  *  TT 


(18) 


where  §2^77  is  linear  in  u^,  v^,  w^ 

is  quadratic  in  these  variables,  etc. 

Substitution  of  the  chosen  forms  for  the  displacements  into  the  ex¬ 
pression  for  the  total  potential  energy  77  and  integration  of  the  resulting 
definite  integrals  is  a  straightforward,  if  tedious  procedure.  The  expres¬ 
sion  for  B^  then  has  the  form: 

^277  =  a^^A^  +  Sa^gAB  +  2a23AC  +  +  SaggBC  +  a^gC^ 


where  the  a^j  will  be  defined  momentarily. 
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The  system  will  be  in  stable  equilibrium  if  Sgvr  >  o,  unstable  if 
8 g7T  <  o.  Hence  the  condition  for  the  onset  of  instability  is 

SgTT  =  o,  (20A) 

or,  in  view  of  the  fact  that  Sgw  is  a  quadratic  form  in  terms  of  A,  B,  C: 


®11 

^12 

^13 

®21 

^22 

^23 

=  o. 

120B) 

^31 

^32 

^33 

Each  of  the  terms  a^j  may  be  written  in  the  form 


ij 


a 


ij 


JL  a 
*  E  ^ij 


(21) 


Hence  the  stability  criterion  is  obtained  from  a  solution  of  Equation  20B 
for  p/E.  Upon  expanding  this  determinant  into  a  cubic  equation  in  p/E, 
it  was  observed  that  the  cubic  and  quadratic  terms  could  be  neglected  com- 
peired  with  the  constant  term  and  the  term  lineeir  in  p/E.  The  resulting 
equation,  when  solved  for  Pqj./E  yields: 


Per  °^ll(°^25^  ~  °^22°^55)  *  *  °^22°-15^  ~  ^12P’15°^g5 _ 

E  Q'1i(®-22^33  *  “33^22^“^11“23^25““12^^55  ^12^“l3^23  “23^13^ 

-  2agga^jj/3^g  -  ^22°^13^  •  (22) 


In  arriving  at  the  expressions  for  the  a^j,  the  following  simplifying 
approximations  were  made:  (a)  Quadratic  terms  in  the  pre-buckling  dis¬ 
placements  Uq,  Vq  and  Wq  (which  are  directly  proportional  to  p/E)  were 
neglected;  (b)  Fourth  degree  terms  in  u^,  v^  and  w-j^  were  neglected; 

(c)  In  view  of  the  form  of  the  chosen  displacements 


W]^  =  C 


fl(/S)  8in(eu'  +  b)cos  n  6  *  fg(/3)sin(ar  +  b)cos  n  6 


. 
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each  is  of  the  form: 

(^ij j)  +  If)  *  ^2^^^ III^ ' 


Now,  apart  from  their  values  as  /3  approaches  zero  and  unity,  the  forms  of 
f^  and  f|g  were  left  unspecified.  We  now  require  further  that  f^  and  fg 
be  chosen  such  that  the  cross-product  terms  become  negligible.  The  amount 
of  error  introduced  by  this  approximation  can  be  systematically  evaluated 
and  reduced  by  the  judicious  choice  of  f-|^  and  fg.  A  little  experimentation 
reveals  that  the  forms  chosen  in  Equation  13  are  quite  satisfactory:  i.  e. 
neglecting  the  cross-product  terns  does  not  significantly  affect  the  deter¬ 
mination  of  the  buckling  load.  We  may,  then,  regard  a^j  to  be  of  the 
form: 


where  the  single  and  double  asterisks  indicate  the  contribution  from  the 
first  and  second  portions,  respectively,  of  the  assumed  displacement  forms. 

Omitting  an  unessential  constant  factor  Eh/2(l-v^),  the  non-vanishing 
and  /3j^j*are  given  by: 


^11  “  "  ^ 


^  h^sin^a 

I*  +  - p-  II* 

12  ro^ 


^12  ”  HTT”  ^ 


^  h^sin^a 

III*  +  - 5“  IV* 

6  ro^ 


'IS 


^22 


TT  sin  a  <  III*  + 


^  |7n^  +  sin^a)  IV*  -  cos^aVI^ 

TT  \  o  ^  .  ,  p  h^sin^a  I~  o 

—  <n2  VII*  +  2  1-v  cos"^aX*  +  - ^  \  XIII* 

sa  I  12  ro'^  L_ 


n 


+  2(1  -  vjcos^aXII* 


®23  “  “  S  VII*  + 


12  ro^  “ 


(25) 


2  +  sin2a)XIII*  -  2  cos^aXIV*  r 


V* 
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ttjjjj  =  -  77 cos  a  tan^a(VII*)  + 


12  ro‘ 


1^0  s^a 


(XI*  +  XII*) 


-  2(n2  +  sin^ajXVII*  +  2^  cos^a  XVIII*  +  8(l-i^)n2  XII* 


^  /(n^  +  sin^a)  ,  .  ,  ^ 

XIII*  - ^ - -  +  8(1 -1^)  (n^  -  2n2) 

\  cos'^a  / 


,*3  .  -  ..A  (^)  V. 


*  77tana  iro\ 
^22  = 


/ro\ 

I'hJ 


'23 


n77 


4co6a 


m 


.-12(l-v)lX*  ,  p, 

VIII*  - ^ H -  +  8(l-i^2)XIX* 

tan'^a 


VIII*  -  6(l-2v)cot2a  IX*  +  8(l-i^2)8ina  XIX* 


(25) 


/S33  =  377cota(l  -  2j^)[— J  XV*  *  2co^  i — 1  Zv  ) 


where  I*  through  XIX*  are  polynomials  in  j3,  and  are  tabulated  in 
Appendix  C(a). 

The  components  a^j  and  arising  from  the  second  portion  of  the 

displacement  functions  may  be  evaluated  in  terms  of  the  Sinus  and  Cosinus 
integrals  as  follows; 


“11  “  - 


TT'^COSa 


f—  -  - 

I2  /3 


/Scosaig 


♦♦  h^sin^a 


12  ro‘ 


/fl 


2~  /Scosaig 


** 


+  77^cosaI*  //3  +  TTCOsalip 


♦♦ 


“12 


'13 


n 


-T  /SI 


P  h^sin%  f”  ** 

*  "  6  ro2  ' 


/Ssinal3  ^'77‘^sina 
2  *  2 

♦ 


h^sina  r”,  p  .  p  , 

- rj-  /S(n^  +  sin^a)l 

24  ro^ 


** 

5 


+  77^(1  +  i^) cos^aIjj//S  +  vTT^cos^a/lfi  +  277lg  cos^a 


(26) 
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a 


22 


+  sin%^ 

p  **„ 

ptana 


•  arlg  ( 

sina 


(l-i^) 


3  /I 

TT'^COSa  — - 

i/3  2) 


h^sin^a 
12  ro^ 


+  |^2(  1  -  i' j^cos®a)/cosa  +  2(l  -  vjTr^cosalg  //S 


** 


+  477(1  -  v)cosaIij, 


“23  "  ^ 


-  ▼▼ 

ptamlQ^  + 


h^sin^g 

4(n^  + 

sin2g);Ql** 

12  ro^ 

sing 

cosg 

)/QcotaIg 


'•33 


O  WW  .  .  P  ▼▼ 

4i/77-'^cotalT  8(1  -  vlTT^cotalo  , 

— ^  ^  ‘ 

+  77cotaI^*(2  +  12[l-^'^ 

**  h^sin^af”  TT^cosa  1  1 

Aanasinal,  .  -  -_ 

/3(n^  +  siting)  ^Ig*  772(co8^a  +  8n^(  1  -  i^) l2*cosa 


ain'^acosa 


/Ssin^g 


2vrf^(n^  *  sin^g)!-!  cosg  8{  1  -  v)n^/ScosgIg 

^  - g  -  2' - -  *  - ^ 

psin'^g 


.  "jr 
Bin'^g 


77cosgl7*(n^[9  -  8v]  +  sin^' 


sin^g 


/S 


**  3 

13  =  “ 


(l  +  v)ti^p\2  ~  /s)(cot^) 


i 


s 


/3 


22 


p» 

^  /ro\  V77(1-2i^)  3(  1  -  v)77COt^g  77(1-^®) 

ptang  —  • -  +  - -  + - 

ihy  4  4  2 


+  (1  -  v/2) 


E 


** 


** 


77  77*^ 


(26) 
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where  through  given  in  terms  of  the  Sinus  and  Cosinus 

Integrals,  and  are  tabulated  in  Appendix  C(b). 


RESULTS  AND  DISCUSSION 

Equations  22,  25  and  26  determine  the  value  of  the  critical  pressure. 
These  relations  were  programmed  for  a  digital  computer  and  Pcp/E  was  deter¬ 
mined  for  the  following  set  of  parameters; 

a  =  15°,  30°,  45°,  60°,  75°,  85° 


h  h  sina 
p  ro 


.02,  .01,  .005,  .002,  .001,  .0005,  .00025 
and  .0001 


^  -  .05,  .1,  .2,  .3,  .4,  .5,  .6,  .7,  .8,  .9,  1.0 
and  v  =  . 3  in  all  cases. 


VIhen  these  parameters  are  inserted  into  Equation  22,  Per/®  still 
remains  a  function  of  n-j^  the  number  of  lobes  in  the  circumferential  direc¬ 
tion.  Since  we  are  seeking  the  minimum  buckling  load.  Pgr/®  minimized 
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with  respect  to  n,  subject  to  the  restriction  that  n  is  an  integer  greater 
than  one.  The  results  of  these  computations  are  plotted  in  Figures  3 
through  8. 

Comparison  of  these  curves  with  the  corresponding  results  of  Radkowski 
(Reference  6)  reveals  that  the  present  results  are,  in  general,  lower  esti¬ 
mates  of  Pqj,/E.  The  difference  between  the  results  diminishes  for  small  /3, 
where  the  conical  frustra  behaves  more  like  a  cylinder,  and  where  Niordson 
(Reference  5)  and  Radkowski  (References  6  and  7)  are  essentially  in  agree¬ 
ment.  The  difference  in  results,  then,  reflects  the  /3  dependence  of  the 
buckled  shape.  Furthermore,  our  results  verify  Radkowski 's  (Reference  6 
and  7)  physical  hypothesis  that  when  plotted  versus  /S,  Pgj./®  does  not  in¬ 
crease  significantly  from  its  minimum  value  (attained  at  or  near  /3  =  .6) 
as  /3  tends  toward  unity. 

The  results  obtained  in  this  report  for  complete  cones  (/S  =  l),  when 
compared  with  the  corresponding  results  of  Grigolyuk  (Reference  8)  are  sub¬ 
stantially  lower  for  larger  values  of  h/p.  As  h/p  decreases,  the  difference 
between  the  results  diminishes,  and  for  very  small  values  of  h/p,  (very  thin 
shells)  the  two  solutions  would  apparently  coincide. 

A  comparison  of  our  results  with  those  obtained  by  Seide  (Reference  9) 
is  complicated  because  of  the  difference  in  the  selection  of  parameters  and 
the  sensitivity  of  the  results  when  interpolation  is  used.  For  this  reason 
we  have  computed  Pq^/E  for  the  same  parameters  employed  in  Reference  9  for 
two  values  of  the  angle:  a  =  85°  and  a  =  30°.  The  ratio  of  our  results  to 
Seide 's  results  for  these  parameters  is  given  in  Table  I. 

Inspection  reveals  that  the  results  of  this  report  seem  to  be  in  closer 
agreement  with  the  results  of  Reference  9  than  other  energy  solutions.  How¬ 
ever  a  difference  in  the  dependence  of  the  solutions  upon  h/p  is  noted:  for 
higher  values  of  h/p  our  results  lie  below  those  of  Reference  9  but  as  h/p 
becomes  very  small,  our  results  become  larger. 

This  difference  in  the  independence  of  P^r/E  upon  h/p  might  well  be 
made  the  basis  of  a  ceureful  experimental  investigation.  It  reflects,  par¬ 
tially,  the  influence  of  the  bending  portion  of  the  strain  energy  which  is 
generally  thought  to  have  much  less  significance  than  the  stretching  energy. 

The  procedure  employed  in  this  report,  then,  when  used  with  suitably 
chosen  displacement  functions,  yields  numerical  estimates  of  critical  pres¬ 
sure  which  tend  to  be  lower  than  those  obtained  by  applying  the  Rayleigh- 
Ritz  method  to  the  classical  strain-energy  procedure. 
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TABLE  I 


=  RATIO  OF  RESULTS  OF  PRESENT  INVESTIGATION  TO  THOSE  OF  REFERENCE  4 


0.0804S 


1.0199 

1.116 

1.177 

1.232 

1.372 


.6787 

.7840 

0.8377 

0.8924 

0.9899 


0.6104 

0.7111 

0.7765 

0.8066 

0.9081 


h/p 


0.0038 

0.00191 

0.00127 

0.0009S4 

0.000477 


0.00366 

0.00183 

0.00122 

0.000916 

0.0004S8 


0.00277 

0.00139 

0.000924 

0.000693 

0.000347 


0.00213 

0.00106 

0.000708 

0.000S31 

0.000266 


.9517 

0.00339 

.054 

0.00169 

.121 

0.00113 

.168 

0.000848 

.317 

0.000424 

■Si 

0.00145 

mggm 

0.000725 

0.6363 

0.000483 

0.9100 

0.000362 

1.0291 

0.000181 

0.7186 

0.8336 

0.9056 

0.9546 

1.0871 


.8547 
.9751 
.045 
1.1045 
1.248 


h/p 


0.00107 

0.000535 

0.000357 

0.000268 

0.000134 


0.000732 

0.000366 

0.000244 

0.000183 

0.0000915 


0.000448 

0.000224 

0.000149 

0.000112 

0.000056 


0.000207 

0.0001045 

0.0000690 

0.00005175 

0.00002588 


0.0001092 

0.0000546 

0.0000364 

0.0000273 

0.00001365 


0.0000562 

0.0000281 

0.00001873 

0.00001405 

0.000007025 
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figure  2 
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F I  CURE  5 
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figure  6 


CRITICAL  PRESSURE  FOR  CONICAL  FRUSTRA 
a  =  75° 


APPENDIX  A 


DERIVATION  OF  STRAIN- DISPLACEMENT  RELATIONS 

The  exact  nonlinear  strain-displacement  relations  for  curvilinear 
coordinates  have  been  formulated  and  are  well  known.  See,  for  example, 
Novozhilov^®  page  59,  where  they  are  presented  in  the  form: 


When  the  usual  thin  shell  orthogonal  coordinate  system  (a-j^,  ug,  z) 
is  used,  and  Ug  being  lines  of  principal  curvature  on  the  shell's 
middle  surface,  and  z  perpendicular  to  both  of  these)  we  have: 


®11 

®22 


1  /  1  Bu  1  dA^ 

1  +  z/R^  \'^1  Bcl A^^Ag  ^^2 

1  /  1  Bv  1  BAg 

1  ^  z/Rg  iAg  ^^2  '^l'^2  Boi^ 


®zz  ®zz 


Bw 

Bz 
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1  3v 


'IS 


u  + 


/  1  3u  1  dA2 


(l  +  z/R^)  v^l  ^  ^  ^  ^^^2  ^"^2  ^^2  '^l'^2 


®13  “  ®lz 


®23  ~  ®2z 


3u  1 

+ 


1  Bw  u 


'dz  1  +  z/R^  \Ai  R-j^ 


Bv 


1  / 1  Bw  V  N 

Bz  1  +  z/R2  I A2  Bag  ^2/ 


Bv  1  /  1  Bw  V  \ 

Sco-.  - - +  -  ( — - 1 

Bz  1  +  z/Rg  \Ag  Ba2  Rgy 


Bu 

Scon  =  - - 


1  Bw  u 


2  “  Bz  “  1  +  z/B^  [Ai  Ba^^  '  R, 


2^3  = 


1  Bv 


BAi 


1  +  ^^^1  \^1  ^l'^2  ^^2 


u 


1  Bu 


BAc 


(A-2) 


where:  components  of  the  rotation  vector  oj,  ^1^2 

the  principal  radii  of  curvature  in  the  a-j^,  a.^  directions,  u,  v,  w  are 
the  displacements  in  the  a^,  ag  and  z  directions,  and  A-j^Ag  are  related 
to  the  Lame  coefficients  of  the  transformation  from  rectangular  to  curvi¬ 
linear  coordinates  by: 

=  Aj  (1  +  z/R-j^)  Hg  —  Ag  (1  +  Z/R2)  • 


By  employing  the  Kirchhoff  hypothesis  (Reference  10)  that  "normals 
to  the  undeformed  middle  surface  go  over  into  normals  to  the  deformed 
middle  surface,  and  undergo  no  extension,  "  and  by  assuming  that  the  dis¬ 
placements  vary  linearly  in  z: 


/s 

U  =  U  + 


z  <p 


V 


A 

V 


+  Z  0 


A 

w  ■=  w  + 


(A-3) 
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it  is  possible  to  express  the  strain-displacement  relations  entirely  in 
terms  of  the  middle  surface  displacements  (u,  v,  w). 

Novozhilov  (Reference  lO)  finds,  for  example,  that; 


where 
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K  (a-4) 


>  (A-5) 


1  Bw 


Substitution  of  these  relations  into  Equation  A-1  yields  what 
Novozhilov  terms  a  "strong  bending  theory”  -  viz. ,  one  in  which  no  assump¬ 
tions  have  been  made  on  the  magnitude  of  deformation  other  than  that  it 
is  one  of  small  strains.  This  formulation,  however,  is  not  only  unwieldy  - 
it  is  unnecessarily  general  for  most  thin  shell  purposes. 
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Most  authors  simplify  the  above  equations  by  using  an  order  of  magni¬ 
tude  argument  in  comparing  the  various  terms  involved  in  the  strain  ex¬ 
pressions.  However,  this  procedure  fails  to  give  proper  attention  to  the 
fact  that  buckling  is  a  geometrically  nonlinear,  if  physically  linear, 
phenomenon.  By  physically  linear  we  mean  that  elongations  and  shear  are 
much  less  than  unity. (small  strains)  and  do  not  exceed  the  proportional 
limit.  This  justifies  the  use  of  a  linear  stress-strain  law  (Hooke's  law). 
By  geometrically  nonlinear  we  mean  that  deformations  may  be  quite  large 
(despite  small  strains)  due  to  the  existence  of  angles  of  rotation  so 
large  that  they  may  not  be  neglected  in  the  determination  of  the  strains. 

Novozhilov  shows  that  when  elongations  and  shear  are  small  compared 
to  unity  and  rotations  are  less  than  unity  but  large  compared  to  elongation 
and  shear,  then  it  is  necessary  to  retain  only  those  terms  of  the  order  of 
the  squares  of  the  rotations  in  the  nonlinear  portion  of  the  strain 
expressions. 


For  thin  shells  [(l  +  z/R-^)  ■=  1,  ( 1  +  z/Rg)  ■=  l] ,  and  assuming  validity 
of  Kirchhoff's  hypothesis  (^22  “  ^ Iz  “  ^2z  “  employ  the  following 
approximate  set  of  strain  displacement  relations.  We  write 


^11  =  ®11  ^  (*^2^  *  "3^) 
^22  =  ®22  ^  ("1^  *  '^3^) 


\  (A-6) 


®12  °  ®12  ■  "l"2  • 


Introducing  Equations  A- 2,  A-3  and  A-4  into  these,  we  have  as  the 
basis  for  an  "intermediate"  bending  theory  of  thin  shells  (einalagous  in 
its  degree  of  accuracy  to  the  nonlinear  VonKarman  plate  equations): 
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where: 


A 

1  bG 

A 

V 

BA]^ 

A 

w 

4.  - 

®11 

AiAg 

Bag 

Hi 

®22 

1  3v 

A 

u 

BA2 

A 

w 

4 

^2  ^2 

A1A2 

Btt]^ 

H2 

®12 

1  3v 

A 

U 

BAq^ 

3ci  ^ 

A1A2 

Ba2 

A 

1  Bu 

A 

V 

BA2 

®21 

A2  ^2 

A1A2 

Bai 

A 

1  3w 

A 

u 

®13 

3cl 

"  % 

A 

1  Bw 

A 

V 

®23 

A^^  "Bcig 

■  _(//+  — 
Ai 

Bw 

A 

V 

2^2 

1 

=  0  - 

Al 

Bw 

Ba 

A 

u 

■% 

2^3 

1  3v 

A 

U 

BA^ 

1  Bu 

V  BAg 

A]^  Ba^ 

A1A2 

Ba2 

A2  ^2 

-^1-^2  Ba^^ 

If 

1  B0 

'P 

4 

ba^ 

+  — . 

*^11 

A^ 

A1A2 

Ba2 

Hi 

If 

1  B0 

4> 

BAj^ 

*^12 

Aq^  3ci 

A1A2 

Ba2 

^^13 

=  ^  ^ 
A^  3ct  ^ 

<P 

'  Hi 

(A-8) 
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^22 


*^21 


*^23 


1  Bt/* 
^2  ^^2 

1 

^2  ^2 

i.  11 

^2  ^°-2 


4>  BAg 

A1A2 

^  BAg 

A1A2 

'/'/R2 


^  :^:7r2 


and,  finally: 


4> 

X 


A 

“  ®23 

®11 


(1  + 
(1  + 

®22 


A  \  A  A 

®22)  ^  ®23®12  “  “  ®13 

A  \  A  A  A 

®ll)  ^  ®13®21  ^  “  ®23 

®11®22  ~  ®12®21  “  ®11  ^  ®22- 


> 


We  observe  that  to  be 
\fj  and  X  since  we  have 
also  neglect  terms  of 


consistent  we  neglect  terms  of  the  order 
already  neglected  terms  of  this  order  in 
the  order  z^T7jj. 


{A-8) 


{A-9) 


4  2  in  0, 

i.  We 
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APPENDIX  B 


DETERMINATION  OF  CHANGE  IN  POTENTIAL  OF  EXTERNAL  FORCES 

For  the  undeformed  cone,  a  point  on  the  middle  surface  of  the  shell 
is  given  by  the  coordinates  (^,^).  These  may  be  related  to  rectangular 
coordinates  by  means  of  the  equations: 

=  (ro  cosa)coB  6 
=  (ro  -  ^  cos  a)  sin  6 

(B-1) 


y 


FIGURE  (B-1) 


For  the  deformed  cone,  a  point  having  original  coordinates  ^,6  has 
final  rectangular  coordinates  X,  Y,  Z  given  by: 

X  =  (r  +  w)cos  6-  V  sin 6* 

Y  =  (r  +  w)sin^+  v  cos6  >  (B-2) 

Z  =  (f  +  u)sin  a . 
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From  these  we  obtain  Gauss'  fundamental  coefficients  E,  F  and  G  by  means 
of  the  relations: 


E  =  +  2^2 

=  {(l  +  u^l^sin^a  +  (w^  -  cosa)^  + 

F  =  (B-3) 

=  (l  +  u^)ugSin^a  +  v^(vg  +  w  +  r)  +  (w^  -  v) (w^  -  cosa) 

G  .  X/  .  y/  .  Z/ 

=  (u^sina)^  +  (r  +  w  +  +  (w^-v)^. 

An  element  of  area  on  the  (deformed)  curved  surface  is  given  by 

dS  =  v/eG-  f2  d^d^  s  Ddfd^,  (B-4) 

hence  the  volume  enclosed  by  the  deformed  middle  surface  is: 

fffdv  =  //Zn^dS  =  ffZn^DdidO  +  J/ZdS  (B-5) 


where  n^  is  the  z  component  of  the  unit  normal  vector  to  the  deformed 
middle  surface. 

But  from  differential  geometry  (Reference  3), 


Dn 


z 


X,Y^  -  Y,X^ 

v^(w^-v)  -  (r  +  +  r  +  w). 


(B-6) 


Now  the  volume  enclosed  by  the  undeformed  conical  shell  is  given  by: 

1  sina  /  , - ) 

Vi  =  -  |Ai  +  Ag  +  v/A^Ag  j  ( 
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where  A-j^  and  Ag  are  the  area  of  the  top  and  bottom  cross  section  of  the 
cone. 


Hence,  the  change  in  volume  during  deformation  is: 

AV  =  JJsina(^  +  u){v^(w^-v)  -  (r  +  w)^(v^  +  w  +  r)}d^d0 

+  sina|lA-]^  +  e(A-]^A2)j-  — +  Ag  +  v/A]^  +  Agj 


\  (B-8) 


where  e  is  the  average  extension  of  the  ends  of  the  cone.  The  work  done 
by  the  uniform  external  pressure  is  then  given  by  p'AV. 
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APPENDIX  C 


(a)  TABULATION  OF  I*  THROUGH  XIX** 

I*  =  /0^(2.  19846)  -^3(10.2868)  +  ^^{20.  IZIO)  -/S(21.6493) 

+  12.3367  -/3"  1(4.  93479) 

II*  =  ^2(4,60870)  -/3(ii.6848)  +  7.40219  -/fl"  1(4. 93479) 

/ 

III*  =  /0^(-.  161731)  +  /33(.  709972)  +  leigi)  +  ^(.942477) 

+  (-.471239) 

IV*  =  /S2(-.  682751)  +  /3(  1.570795)  +  (-1.570795) 

V*  =  y35(-l.  159947)  +  /S^( +6.  903993)  +/fl^(- 20.  51448)  +  y32(30.  16964) 
+/3(-24.  7400)  +  (9.89600) 

VI*  =  /32(4.801058)  +  /S(-16. 82460)  +  (21.47553)  +  (j)(-9. 30186) 
‘^)  (+65093) 

VII*  =  /3^(. 087005)  +  /3^(-.  424009)  +  /8^(.75)  +  /fl(-.5) 

VIII*  =  /3^(-.  250022)  +  /3^(1.  50159)  +  /fl^(-3.  62159)  +  /fl^(4.  24) 

+  /3(-2.  12) 

IX*  =  ^^(.0285192)  +  /34(-.  174009)  +  424008)  +  /3^(-.5) 

+  /S(.25) 

X*  =^(.910570)  +  /3^(-5.3610)  +  /S2(  13. 08704)  +  /3(-16.  9493) 

+  12.33699  +  (1//3)  (-4.  93479) 
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XI*  =  /82(36.8676)  +  /8(-162.  955)  +  (292.0113)  +  ( l/;Q)  (-275. 848) 

+  (1//82)  (121.  761)  +  (l//8®)  (-48.  7044) 

XII*  =  /82(i.  608703)  +  /3(-5.  68480)  +  (7.40219)  +  ( l//8)  (-4.  93479) 
XIII*  =  /82(.25)  +  (;3)(-.5) 

XIV*  =  /82(-2.  059795)  +  >3(8.  314152)  +  (-12.38122)  +  ( l//8)  (8.  38915) 

XV*  =  /35(-i.  79703)  +  /S4(12.  49638)  +  /S®(-37.  3091)  ♦  /8®(61.  96213) 

♦  y3(-61.  82351)  +  (37.01096)  +  ( 1//8)  (-12.  33699) 

XVI*  =  y35(-.  625055)  + /S^( 3.  75397)  +  /S®(- 9.  053968)  +  fi^{l0.60) 

+  y3(-5.30) 

XVII*  =  y32(_.  157609)  +>3(1.05544)  +  (-2.220657)  +  ( l//fl)  ( 1.48044) 
XVIII*  =  /S2(4,  93479)  +  ^(-19.73919)  +  (29.60876)  +  ( l//fl)  (-19.  73917) 

XIX*  =  /35(-.  171116)  +  >34(1.044054)  +  >83(-2. 54405)  +  fi^(5.0) 

+>3(-l.  5). 


-  J 


(b)  TABULATION  OF  I^^*  THROUGH  In** 
■n  sin^  dx  d 


o  1  -  ax 


log(l-a7r)  -  cos2d  Ci  (Y)  +  8in2d  Si  (Y)J 


2b 

2d 


T  ** 
■‘■2 


,7r  COS^X  dx  ,  ,  .  .  . 

— - -  =  -  d  log  ( 1  -  an)  -  1-^** 


o  1  -  ax 


=  J 


■n  sin2x  dx 
o  1  -  ax 


=  -  d  [sin2d  Ci  (Y)  -  cos  2d  Si  (Y)] 


2b 

2d 


Ia**  =  J 


tr  COS 


dx 


o  ( 1  -  ax ) 3 


_d 

2 


( 1  -  an)  2 


-  1 


T  ** 
■‘■6 
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Is"  -  / 


77  sinSx  dx 
o  ( 1  -  ax) 3 


1 

_d  '  b_  ~ 


2  T  ** 


Bd'^  I3 


■r  **  r'"  sin^x  dx  o 

^6  =  f  M  =  -  d^logd-^)  - 


o  ( 1  -  ax) 3 

p77  sinSx  dx 
o  ( 1  -  ax)  2 


,77  sinSx  dx  o 

^7  “  -^n  7TT~r2  “  log  ( 1  -  arr)  +_  4d  ** 


^77  sin^  dx 

■f  T. - ^  ■=  -  ^13 

o  ( 1  -  ax)  2  o 


rft  cos^  dx  /  v\ 

4  71^2  ■  *t) 


T  ^** 
-‘■10 


,77  xsin^dx  o 

J  7- - ^  =  -  dl,**  -  d2  I3** 

o  (1  -  ax)2  ■>■  o 


,  **  r’'  x^sin^  dx  p  a-a.  9 

'11  ■  4  (1  -  ax)2  °  "  ^3"  *7<i" 


where : 


1  TT 

d  =  —  =  —  b=d-7T 
a  yS 


,  ^  2b 
Ci  Y) 

■  X 

cos  y 

dy 

2d 

2d 

y 

...  2b 
Si  Y 

«  : 

III 

sin  y 

dy  . 

2d 

2d 

y 
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